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Abstract 

The heavy quark potential and particularly the one proposed by Richardson 
to incorporate both asymptotic freedom and linear confinement is analyzed in 
terms of a generalized Borel Transform recently proposed. We were able to 
obtain, in the range of physical interest, an approximate analytical expression 
for the potential in coordinate space valid even for intermediate distances. The 
deviation between our approximate potential and the numerical evaluation of 
the Richardson's one is much smaller than A of QCD. The cc and bb quarkonia 
energy levels agree reasonably well with experimental data for c and b masses in 
good agreement with the values obtained from experiments. 



Classification: 12.38 Aw - 12.38 Lg 

Keywords: quark-antiquak potential; Borel transform; analytic properties. 



Partially supported by CONICET and ANPCyT-Argentina. 



1 



Among the different proposals for describing quark-antiquark interactions, the Richard- 
son's potential due to its is simple structure, has been the subject of continuous 
interest 0. This potential requieres, in the author words, the minimal number of pa- 
rameters. In fact, the only parameter entering the potential is the QCD related scale 
A. This potential, designed in order to present both asymptotic freedom and linear 
quark confinement includes the single dressed gluon exchange amplitude, namely 



v(^)-l"-jf^ (1) 



Asymptotic freedom is present as soon as one adopts for as{q^) the effective run- 
ning coupling constant provided by the renormalization group. Quark confinement is 
imposed by requiring that for q small V{q) behaves as the inverse four power of q that 
guarantees a linear behavior in r. Then, the Richardson's potential reads 

V^(q) = (2) 

where Cp is a. group coefficient. 

It is clear that the explicit calculation of the QCD coupling constant in position 
space is crucial when the Richardson's potential is to be applied in a concrete case. 
This is because the Fourier transform of @ provides the configuration space expression 

V (r) = — a{l/r) = — a(l/r) - A r (3) 

r Po Po 

Here 

a(l/r) = l-4/(r) (4) 

f(r) = r ^ e-g^^ 

J I q [ln(g2 - l)]2 + ^ ' 

This expression was only computed numerically. There exist some analytical results 
corresponding to some asymptotic conditions. For example, for Ar <^ 1, the Richard- 
son's potential was shown to behave softer than the Coulomb interaction [Q, namely 

^""(^^ ^ Arln(Ar) 

and for Ar ^ 1 provides linear confinement. 

Our main point in this paper is the calculation of the strong coupling constant in 
position space starting from the last integral representation eq.(||). In so doing, we 
provide either an input for the Richardson's potential in configuration space or to any 
other alternative proposal for the quark potential including the original QCD running 
coupling constant or any alternative expression ||^. To this end, we fully analyze the 
analytic structure of the integral in the Borel plane [^]. Then, we are able to obtain 
the potential behaviour as a function of r, including intermediate distances. 
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Any perturbative analysis starts from the general relationship between V{q) and 
V{r) that ends with the corresponding relation between the couplings as{q) and a(l/r). 
Notice that in the perturbative calculation, a(l/r) coincides with a(l/r). Moreover, it 
has been shown [0 that one can write, for any static potential. 



a(l/r) =J2fn 



d 

-(3{a ' 



da. 



a,{q = K/r) (7) 



where /„ are known constants, k = exp{'j^) a constant and l3{as) = fi^ das{q)/dfi'^ . In 
the case of the Richardson's potential, this series is asymptotically (g^ ^ 1) factorial 
divergent and its Borel sum does not exist, namely 

a(l/r) ~ as{q = n/r) J^fn [/5o"^(9 = '^A)]" ^' 

n 

Certainly, this expression provides sensible results only for very small distances because 
at increasing distances the non-perturbative contributions start to be important. 

It is worth mentioning that the analytic behaviour of a{l/r) has been studied 
by summing the divergent asymptotic series by using the standard Borel formalism. 
Clearly, being the expression no Borel summable, the Borel transform B{s) has singu- 
larities for different values Sk on the integration path of 

POO 

a{l/r) = / exp[—s/a{q = K/r)]B{s)ds 
Jo 

Consequently, it can be defined only in principal value, showing ambiguities coming 
from the exponential in this integral. This approach implies that the non-perturbative 
contribution is considered of the same order of magnitude as the ambiguity inherent to 
the method [^]. An additional problem coming from the use of a perturbative is the 
presence of the Landau pole, conditioning the validity of any amplitude representing 
any physical observable to a finite range of energy. In this respect, there is an alter- 
native proposal 1^ that starts from a modification of the definition that avoids the 
Landau pole but retains the standard properties. Nevertheless, this change implies a 
modification in the linear confinement behaviour loosing the standard connection with 
the string tension. 

All the previous mentioned problems can be avoided by using the generalized Borel 
Transform (GBT) that was introduced in Ref. |TU|. This version of the Borel transform 



was originally defined on a finite lattice but it can be readily adapted to the continuum, 
preserving all of its characteristics. 

The main vantage of this proposal comes from the fact that its analytic properties 
have no ambiguities. Moreover, it allows to perform computations in terms of a real and 
positive arbitrary parameter A avoiding the implementation of perturbative expansions. 
The approach generally ends with non-perturbative calculations of the saddle point 
type, when A takes large arbitrary values. This is possible because the generalization 
implies the definition of a valid Borel transformation for each value of the parameter 
A. Then, the better adapted value for each particular problem can be chosen. In other 
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words, when using the GBT, a function of Bx{s), a whole class of transformations is 
performed. It is found that, as it should be, the results do not depend on A. For the 
particular case of the present paper, this can be summarize as 

a{l/r) =T^^[Tx{a{l/r))] where Tx{a{l/r)) = Bx{s) for < r < oo 

We start by presenting our previous generalization [|10| of the Borel transform of a 
function /(r), namely 



Bxis) 



exp[s / ri{r)] 



Xri{r] 



+ 1 



-As 



/(r)d(lMr)) ; Re{s) < 



(8) 



Among the properties of this definition we want to notice that it is valid for any analytic 
function ri{r) in the interval < r < cxo that allows to define 



1 



uxir 



- A In 



1 



A ri{r] 



+ 1 



(9) 



being monotonically increasing in the same interval if ri{r) is. Consequently, the integral 
transform can be written as 



Bxis) 



exp[s u] f[rx{u)] {1 + \ri[rx{u)]} du ; i?e(s) < 



(10) 



where rx{u) is the inverse coming from the change of variables. From the last expression 
it is clear that 



Bxis) 



POO 

/ exp{su) Lx\r{u)\du ; Re{s) < 
Jo 



(11) 



is the Laplace transform of the function Lx[r{u)] implicitly defined. That definition 
implies that r]{r) gives rise to an analytic transformation in the negative Borel half- 
plane, such that its extension to the other half-plane is also analytic with a cut on the 
real positive axis. From this observation it is clear that /(r) can be unambiguously 
expressed in terms of the inverse Laplace transform integrated on the above mentioned 
cut (for details see Ref. |llO| ) 



fir) 



A rjir) 



exp[—s/riir)] 



A rjir] 



+ 1 



Xs 



ABxis)ds 



(12) 



As it was said before, the parameter A can take any real positive non zero value 
generating a continuous family of transformations. A large value of A could be useful 
because in this case asymptotic techniques can be used in the calculations. 

From this point on, a series of almost trivial calculations follows. Let us only 
indicate the most important steps. Using the ansatz l/?7(r) = A [exp(Ar/A) — 1] 
which is well defined for < r < oo, the function uxir) results 



Uxir) = A [expiAr/X) — 1] — A r 
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Consequently, one can write 



with 



dqH{q) 
H{q) 



exp{sXv) {v + l)-^^'+'''^dv 



(13) 



(14) 



[In(g2-l)]" + 7r2 q 

where the change of variable l/[Xri{r)] = v was introduced. Notice that the last inte- 
gral, for Re{s) < 0, represents the confluent hypergeometric function [TT[| G{1, —As — 
Xq + 2, —s X). Consequently, in this region Bx{s) is an analytic function and when 
an analytic continuation to the positive half-plane of s is performed, a cut appears. 
Introducing the discontinuity of the G function, one gets 



AB,is) = 2nXexp{-Xs) 

We can now transform back to obtain /(r) from eq. (|l2|) 

1 



(15) 



/ (r) = - A, (r) 

TT 



exp [G {w, t, r. A)] dw dt 



(16) 



where 
and 



Ax (r) = 1 - exp [-Ar/A] 



G{w,t,r,X) = — Xv (t) ux{r) +t — nw — In 



w 



2 In (g (w)) 



-\n{r[X{Xv{t) + q{w))]} + Xq{w) In [X^ v {t)] - X^ v (t) 
+ [X^v{t) - l) In [X^v{t) 



[17) 



with 



q[w) 



l + e 



1/2 



vit) 



The next step is to look for the asymptotic contribution in A of the double integral 
in eg. ([TBI) . To this end one can use the steepest descent technique in the combined 
variables {t,w). Consequently one first computes the saddle points to^r) and Wo{r) 
and then one checks the positivity condition in particular when the discriminant 
D{tQ, Wo) of the second derivatives of G at this point is positive. In so doing one obtains 



/ (r) ~ 2 A' Ax (r) exp [G {wo (r) , to (r) , r. A)] 



d^G d^G 



dw"^ dt"^ \dwdt J 



-1/2 



(18) 



the saddle point being 

<lo{r) 



to = In 



Fiqoir)), 



to{r) 



Wo 



In 



9o W 



Wo (r) 



(19) 
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where qo{r) is the solution of the imphcit equation coming from the extremes of the 
function G, namely 



A2 



with 



Fiqo) 



2qi\n{qi-l] 



(20) 



(21) 



go [go'-i] r [In(go'-i)]' 

Notice that -F(go) should be positive and consequently go < 2.130156. On the other 
hand, from eq.(|l9l), go > 1. In fact, moving go between these values, the variable r 
covers all the positive real axis in a biunivocal way. Moreover, the condition -F(go) 7^ 
implies that r = is excluded from the analysis. This is clearly not a drawback of the 
method. 

Going now to the expression (|l^) of f{r), one finally finds, in the saddle point 
approximation (A ^ 00) 



fAp{r)=f{qo{r)) 



where 



e-'/'V2j{F{qo) + l^ 



2j7rD{qo 



[go 



1 3/2 



go^ 



exp [-go F (gp)] 

^ln(go^-i; ^ 



(22) 



^(go) 



2g 



^ Q'oF{qo) 



1 - 



+ - 



^in(go^-i) 



^ In (gg - 1) 



+ 1 



(23) 



TT (gg -1)F ( go) 

2go 



1 2 



Finally, the approximated expression for the potential is 



VAp{r) = - 



2 7rCf 



(24) 



In order to test the precision provided by the GBT, we have compared the behavior 
in the coordinate space of our approximated analytical formula (24) with the numerical 
integration of the exact expression (3). This comparison, within the range of physical 
interest 0.1 /m < r < 1 fm (see for example |T^, 0), has been performed for the 
Richardson's potential free parameter A = 0.398 GeV which he found to be the best to 
fit the charmonium and botomonium lowest bound states. The corresponding results 
are exhibited in Table I, being the relative error 5% at most. 

We have also evaluated, the charmonium and bottomoniun spectra following the 
Richardson's algorithm. For this purpose, with the help of (24), we have adjusted rric 
and A to reproduce the best measured charmonium energy levels ip{lS) and 4'{2S) [H . 
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The resulting charm quark mass is rric = 1.363GeV^ and the scale size A = OASAGeV (to 
be compared with = lA91GeV and A = 0.398Gel^ obtained by Richardson). In the 
bottomonium case we kept A = 0A34GeV assuming, like Richardson, that the scale size 
should be a fundamental scale of the theory independent from the quarkonium system. 
The bottom mass reproducing the best measured bottomonium energy state T(25') 



14| results in mf, = 47806*6^^ . Notice that Richardson has fitted the botomonium 
fundamental level T(15') obtaining = 4888G'el^. The results for both spectra and 
the corresponding experimental values are presented in Tables II and III respectively. 

In Summary, we have obtained an analytic expression for the quark potential valid 
for any value of r. Our results, in the leading order of the saddle point approximation, 
give rise to absolute values of the uncertainties always much smaller than A of QCD 
(~ OAGeV). It is worth mentioning that any perturbation based calculation obtained 
by means of the standard Borel transform, ends with uncertainties at least of the order 
of A 1]^. This result shows the ability and potentiality of the previously introduced 
generalized Borel Transform. 

As regards the quarkonia energy levels, our spectra prediction results reasonably 
well when compared with the spectra obtained by Richardson. The deviations of our 
results with respect to the corresponding experimental values, are certainly smaller 
than the scale size A (see Tables II and III). Let us remark that our adjustment leads 
to quark masses, when compared with those of Richardson, in best agreement with the 
values reported in the particle data table [0]. The resulting A parameter falls in the 
range of values expected from perturbative QCD. 

There exist in the literature other simple potentials in the coordinate space P ||T3 



that can compete with our approximation and that give similar results for the spectra. 
However these potentials are strictly phenomenological, not having a direct connection 
with pertubatibe QCD and in general, they have more free parameters to be adjusted. 

Finally, our QQ approximate potential, due to computing simplicities, is clearly 
useful in any further calculations of interesting physical quantities |16[ . 
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Table I: Comparison between potentials for A = 0.398GeV 





' Rich [^<~- y J 


Va \GeV] 


\()V\ 


7.0 


0.899941 


0.899923 


0.000018 


5.0 


0.554276 


0.554244 


0.000032 


3.0 


0.149355 


0.149637 


0.000281 


1.0 


-0.571347 


-0.554343 


0.017004 


0.5 


-1.114402 


-1.062477 


0.051926 


0.2 


-2.208255 


-2.284687 


0.086730 



Table 


II: System cc 






/ 


Ee.p [GeV] 


Emch [GeV] 


Eap [GeV] 





3.096 


3.096 


3.096 





3.684 


3.684 


3.684 





4.040 (4160)* 


4.096 


4.127 





4.415 


4.440 


4.506 


1 


3522 


3.514 


3.494 


2 


3770 


3.799 


3.786 



Our parameters are [A = 0.434GeV, rric = 1.363GeV] and the corresponding Richard- 
son's ones [A = 0.398GeV,mc = lA91GeV]. 
* See l]!^ for an alternative level assignment. 

Table III: System bb 



I 


Ee.p [GeV] 


Emch [GeV] 


Eap [GeV] 





9.460 


9.460 


9.522 





10.023 


10.016 


10.023 





10.355 


10.343 


10.351 





10.580 


10.607 


10.620 


1 


9.900 


9.896 


9.910 


1 


10.260 


10.249 


10.252 



Our parameters used are [A = 0. 4346*6 V, = A.780GeV] and the corresponding Richard- 
son'sq ones [A = 0.398GeV",m6 = 4.888GeV]. 
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